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We derive stability conditions of Asymmetric Nuclear Mat- 
ter (ANM) and discuss the relation to mechanical and chem- 
ical instabilities of general two-component systems. We show 
that the chemical instability may appear as an instability 
of the system against isoscalar-like rather than isovector-like 
fluctuations if the interaction between the two constituent 
species has an attractive character as in the case of ANM. 
This leads to a new kind of liquid-gas phase transition, of in- 
terest for fragmentation experiments with radioactive beams. 

PACS numbers: 21.65. -hf, 25.70.Pq, 71.10.Ay 

A binary system manifests a richer thermodynamical 
behaviour as a consequence of a new parameter, the con- 
centration, which is required for a complete description 
of its states. The phase transitions are more complex be- 
cause they have to accommodate one more conservation 
law. A thermodynamical state can be unstable not only 
mechanically but also chemically. 

Nuclear matter belongs to this class of systems. The 
process of multifragmentation following the collision of 
heavy nuclei at medium energies is expected to show fea- 
tures analogous to usual liquid-gas phase transitions of 
water ||^-^. In particular, the spinodal decomposition 
appears to be an important mechanism leading to phase 
separation in symmetric nuclear matter (SNM) ^^, 
as confirmed in recent experiments [J7|. The relevance 
of instabilities in AA^M-fragmentation has been already 
discussed few years ago |^-[l(]|| . It was shown that in these 
systems a kind of diffusive (or chemical) spinodal rather 
than the mechanical spinodal is significant. However the 
detailed nature of fluctuations responsible for such insta- 
bility has not been clarified yet. This is a quite important 
aspect since, as we will see, it provides reliable informa- 
tion on the character of the interaction in the medium. 
We remark that the effect is leading to the "isospin dis- 
tillation" that can be experimentally observed, see the 
recent data |ll|. 

In one-component systems the mechanical instability 
is related to instability against density fluctuations as 
the result of the strong attraction between constituents. 
In symmetric binary systems, like SNM, one encoun- 
ters two kinds of density fluctuations: i) isoscalar, when 
the densities of the two components oscillate in phase 
with equal amplitude, ii) isovector when the two den- 
sities fluctuate still with equal amplitude but out of 
phase. Mechanical instability is associated with insta- 



bility against isoscalar fluctuations leading to cluster for- 
mation while chemical instability is related to instability 
against isovector fluctuations, leading to species separa- 
tion. We will show that in asymmetric binary systems, 
as in ANM, this direct correspondence between the na- 
ture of fluctuations and the thermodynamical instability 
is lost. 

An appropriate framework for this investigation is pro- 
vided by the Fermi liquid theory which has been already 
applied to the study of instabilities in symmetric binary 
systems as SNM (the two components being protons and 
neutrons) and the liquid ■^iJe (spin-up and spin-down 
components) |12|Jli^ . 

We first investigate the thermodynamical stability of 
ANAI at T = extending to the asymmetric case the 
formalism introduced in [[l2[ . The distribution functions 
for protons and neutrons are: 

/f (6^ = 9(^^,-6^, q = n,p (1) 

where fj,q are the corresponding chemical potentials. The 
nucleon interaction is characterized by the Landau pa- 
rameters: 



pqiq2 = Nq,V^ 
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where Ti. is the energy density, V is the volume and Nq is 
the single-particle level density at the Fermi energy. At 
T = this reduces to Nq{0) = mpF,q/{Tr'^ti^), were pp.q 
is the Fermi momentum of the g-component. Thermo- 
dynamical stability for T = requires the ground state 
energy to be an absolute minimum for the undistorted 
distribution functions, such that the relation: 



SH - ^IpSpp - finSpn > 



(3) 



is satisfied when we deform proton and neutron Fermi 
seas. The distorted distribution functions can be written 
as /g(p) — 6(6F,g(^)— 6^) , where eF,q{d) is a direction de- 
pendent Fermi energy characterizing the distortion. We 
will follow the usual multipole expansion for the variation 

[SeF,,{0)] = ^F,,(e) 'fiq^Y. ^[Pii^osO) (4) 

and for the Landau parameters 
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We limit ourselves to monopolar deformations consider- 
ing here momentum independent interactions such that 
^/=o^ are the only non-zero Landau parameters [[l6[ . 
Then, up to second order in the variations, Eq.(y) be- 



rhei 
comes 



SH - fipSpp - ^nSpn = ■:^ia'^p + bl^n + CVpVn) (6) 

where Vn,p = i^n.p ^^'^ 

a^Np{0)il + FPn ; 5 = iV„(0)(l + i^^"") ; 

c = 7V„(0)i^P" + NpiO)F:;P = 27V„(0)Fr. (7) 

We diagonalize the r.h.s. of Eq.(g) by introducing the 
following transformation: 



u — cos [3 Up + sin(3 v^n 
V — ~sinj3 Vp + cos (3 Vn, 

where the mixing angle < /3 < 7r/2 is defined by 



(8) 
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c _ JV„(0)fo^"-HiV^(0)fo"P 

. - b Np{Q){l + F7) - 7V„(0)(f + Fg"") ■ 



Then Eq.(p[) takes the form 

5T-i — fipSpp — llnSpn = Xu^ + Yv^ 

where 

1 
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X = -( a + b + sign{c)^/{a - 6)2 + c? ) 



and 



A^p(0) + iV„(0) 



1 



( 1 + FL ) 



(11) 



F = -( a + 6 - sign{c)^/{a - b)^ + c^ ) 
_ JVp(0) + JV„(0) . , ™ . 

= ^ ( 1 + ^Og )• 



(12) 



Thanks to the rotation Eq.(||) we separate the total 
variation Eq.(ph into two independent contributions, the 
"normal" modes, characterized by the "mixing angle" /3, 
which depends on the density of states and the details of 
the interaction. 

In the symmetric case, Np — Nn = -/V, Fq'" = Fq^ and 
pnp ^ ppn^ ^^j^ reduces to 



STL — flpSpp — flnSpn 



NiO) 
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where F^ = F^^' + Fq^ and F^ = F^^^-F^^ are symmet- 
ric and antisymmetric (or isoscalar and isovector) Landau 
parameters, and we recover the usual Pomeranchuk sta- 
bility conditions for pure isoscalar/isovector fluctuations 

In the general case we interpret u- and w-variations as 
new independent isoscalar-VikB and isovector-\\ke fluctu- 
ations appropriate for asymmetric systems. The proton 
and neutron densities will fluctuate in phase for isoscalar- 
like variations and out of phase for isovector-like varia- 
tions, see Eq.(||). Moreover, FJ^ and Fg^, defined by Eq. 
(|ri|,^, can be considered as generalized symmetric and 
antisymmetric Landau parameters. 

From Eq.(|lO|) we see that thermodynamical stability 
requires X > Q and F > 0. Equivalently, the following 
conditions have to be fulfilled: 



1 + FJg > and 1 + F^g > 0, 



(14) 



They represent Pomeranchuk stability conditions ex- 
tended to asymmetric binary systems. 

The new stability conditions, Eq.([l4|), are equivalent to 
mechanical and chemical stability of a thermodynamical 
state, p3, i.e. 



^) >0 and f^) >0 (15) 

dpJr^y \dyJTp 



where P is the pressure and y the proton fraction, as it 
can be proved by observing that ||l8|| : 

XY - 7Vp(0)7V„(0)[(l + F™)(1 + F^P) - F^F^] 
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From Eq.s( [l0| , [l4|) , we are led to define as isoscalar 
instability the case when the state is unstable against 
isoscalar-like fluctuations, i.e. when 1 -I- Fq < (or 
X < Q). Analogously we deal with isovector instability 
when the system is unstable against isovector-like fluctu- 
ations i.e. when 1 -I- Fq^^ < (or 1" < 0). 

One can easily see that in symmetric nuclear mat- 
ter the isoscalar instability, {X < 0,Y > 0), appears 
as mechanical instability and the isovector instability, 
{X > 0,Y < 0), as chemical instability. Indeed now 
t = 1, o = 6, /3 = 7r/4 and so 1 + Fg^ (or X) and 



^ J are proportional (see Eq.(|l7|)). This is not true 

any longer for an asymmetric system. 

To simplify the analysis let us assume for the moment 
that the quantities a and b remain positive in the density 
range we discuss here and so we can study the effect of 
the interaction between the two components, given by c, 
on the instabilities of the mixture. This is indeed the 
case of nuclear matter. 

If c < 0, i.e. for an attractive interaction between the 
two components, from Eq.(^2|) we see that the system is 
stable against isovector-like fluctuations but it becomes 
isoscalar unstable if c < —2\/ab (see Eq.([ll|)). How- 
ever thermodynamically this instability against isoscalar- 
like fluctuations will show up as a chemical instability if 
(— ia — b/t) < c < —2\/ab or as a mechanical instability 
if c < (-ia - b/t) < -2Vab (see Eq.(|l^)). Therefore 
isoscalar-like instabilities are not necessarily equivalent 
to mechanical instabilities and may manifest instead as 
chemical instabilities. 

If c > 0, i.e. when the interaction between the com- 
ponents is repulsive, the thermodynamical state is al- 
ways stable against isoscalar-like fluctuation, but can be 
isovector unstable if c > 2\fah. Since the system is me- 
chanically stable {a,b,c > 0, see Eq. (O)), the isovec- 
tor instability is now associated with chemical instability. 
Such situation will lead to a component separation of the 
liquid mixture. 

Following this line a complete analysis of the instabili- 
ties of any binary system can be performed, in connection 
to signs, strengths and density dependence of the inter- 
actions. 

We show now a quantitative calculation for asymmetric 
nuclear matter which illustrates the previous general dis- 
cussion. We investigate the instabilities of ANM charac- 
terized by a potential energy density of Skyrme-like type. 



'HpotiPn, Pp) 



A {pn + Pp)'^ , B (p„ -f Pp) 
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Po a + Z 

+ (Ci-c2(^r) ^^"~^^^' 

Po Po 



p^' 



(18) 



now the instability can manifest either as chemical or 
mechanical depending on the relative strength of the in- 
teractions in the system, as discussed previously. 




where po = 0.16 fm is the nuclear saturation density. 

The values of the parameters A = —356.8 MeV, B = 
303.9 MeV, a = 1/6, Ci = 125 MeV, C2 = 93.5 MeV 
are adjusted to reproduce the saturation properties of 
symmetric nuclear matter and the symmetry energy co- 
efficient. We also extend the discussion to finite temper- 
atures p5[ . 

We first focus on the low density region, where phase 
transitions of liquid-gas type are expected to be seen 
in fragmentation events. Since a, 6 > and c < 0, in 
the case of the liquid-gas phase transition in asymmet- 
ric nuclear matter we deal only with instability against 
isoscalar-like fluctuations as it happens for symmetric nu- 
clear matter. However, at variance with the latter case, 
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FIG. 1. Spinodalline corresponding to isoscalar-like insta- 
bility of asymmetric nuclear matter (circles) and mechanical 
instability (crosses) for three proton fractions: y — 0.5 (a), 
y = 0.25 ■(b),y = 0.1 (c). 

In Figure 1 the circles represent the spinodal line cor- 
responding to isoscalar-like instability, as defined above, 
for three values of the proton fraction. For y < 0.5 under 
this border one encounters either chemical instability, in 
the region between the two lines, or mechanical instabil- 
ity, under the inner line (crosses). The latter is defined 

by the set of values {p,T) for which (§^) =0. We 

observe that the line defining chemical instability is more 
robust against the variation of the proton fraction in com- 
parison to that defining mechanical instability: reducing 
the proton fraction it becomes energetically less and less 
favoured to break in clusters with the same initial asym- 
metry. 

From Eq.s([l^) and (|l^), we see that the quantity 
( "3lf ) changes the sign passing through zero when we 

cross the spinodal line corresponding to chemical instabil- 
ity and passing through infinity when we cross the line 
associated with mechanical instability. However, from 
Eq.(p^ (see also Fig. 3), we get that at the inner spin- 
odal line, the ratio ( -^ ] / I -^ I oc XY is a finite 

negative number. Therefore we conclude that the singu- 
larities of two quantities cancel by taking the ratio. Note 
also the smooth behaviour of the mixing angle j3 through 
the mechanical instability line. Fig. 2a. 
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FIG. 2. Density dependence of the mixing angle /3, Eq.tel), 
(a) and of the function x (b) for three proton fractions, y — 0.5 
(solid), y = 0.4 (open circles), and y = 0.1 (full circles) at 
T=lMeV. 

Finally let us stress that, for ANM, even if we deal 
with instabilities against isoscalar-like fluctuations, we 
encounter a chemical effect, which gives raise to isospin 
distillation in equilibrium phase transitions |ll| ]. Indeed 
the variation of the asymmetry (7=1 — 2y) is 

^^ = "P^ITTT^^^P - 1] ^ ^p[x - 1] (19) 

11 + lo)I\„ 



strength between the various interactions acting in the 
system. 

In particular, the liquid-gas phase transition in asym- 
metric nuclear matter results from instability against 
isoscalar-like modes rather than isovector-like due to the 
attractive character of the interaction between protons 
and neutrons. This is a qualitatively new effect which is 
leading to the observed "isospin distillation" ||ll|. More 
data are expected from the new Radioactive Beam facili- 
ties opening the possibility of direct studies of the charge 
dependent part of the nuclear interaction far from normal 
conditions. 



where Iq is the initial asymmetry. For y — 0.5 we have 
51 = but for y < 0.5, we find x < 1 (see Fig. 2b). 
Therefore SI < if i^p > and the opposite for i^p < 0. 

To be more general, we have extended our investiga- 
tion also to the high density region. In Fig. 3 we plot 
the density dependence of the generalized Landau pa- 
rameters. We find that for the considered interaction 
the system exhibits another instability at high density, 
around 1.5/m~'^, where the quantity c becomes positive 
and l + F^g <0 (or r < 0). 
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FIG. 3. Density dependence of the generalized Landau pa- 
rameters for two proton fraction, y = 0.4 (a) and y = 0.1 (b) 
(symmetric, solid and antisymmetric, dashed) at T=lMeV. 

Thermodynamically we see from Eq.(|lq) that this is 
again a chemical instability. However now it results from 
isovector-like fiuctuations, in contrast to the low density 
instability. The reason is the change in the character 
of the interaction between the two components. Since 
the interaction becomes repulsive the nuclear phase can 
become unstable against proton-neutron separation. We 
also notice in Fig. 3 that the generalized Landau parame- 
ters display a discontinuity where the quantity c changes 
the sign. Other effective forces, with more repulsive sym- 
metry terms, will not show this high density chemical in- 
stability [|7| , that actually could be of interest for other 
many body systems. 

In conclusion, we have shown that in asymmetric bi- 
nary systems the relevant instability regions are defined 
by the instabilities against isoscalar-like and isovector- 
like fluctuations. The kind of thcrmodynamical instabil- 
ity, chemical or mechanical, will depend on the relative 
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